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ABSTRACT 

Aims. We present the first Doppler imaging study of the two rapidly rotating, single K-giants DP CVn and DI Psc in order to study 
the surface spot configuration and to pinpoint their stellar evolutionary status. 

Methods. Optical spectroscopy and photometry are used to determine the fundamental astrophysical properties. Doppler imaging 
is applied to recover the surface temperature distribution for both stars, while photometric light-curve inversions are carried out for 
studying the long-term changes of the surface activity of DP CVn. Surface differential rotation of DP CVn is estimated from cross- 
correlating the available subsequent Doppler reconstructions separated by roughly one rotation period. 

Results. Both stars appear to have higher than normal lithium abundance, LTE logn of 2.28 (DPCVn) and 2.20 (DIPsc), and are 
supposed to be located at the end of the first Li dredge-up on the RGB. Photometric observations reveal rotational modulation with 
a period of 14.010 d (DPCVn) and 18.066 d (DIPsc). Doppler reconstructions from the available mapping lines well agree in the 
revealed spot patterns, recovering rather low latitude spots for both stars with temperature contrasts of AT « 600 - 800 K below the 
unspotted photospheric background. Spots at higher latitudes are also found but either with less contrast (DP CVn) or with smaller 
extent (DIPsc). A preliminary antisolar-type differential rotation with a = -0.035 is found for DPCVn from cross-correlating the 
subsequent Doppler images. Long-term photometric analysis supports the existence of active longitudes, as well as the differential 
rotation. 

Key words, stars: activity - stars: imaging - stars: late-type - stars: starspots - stars: individual: DPCVn- stars: individual: DIPsc 
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1 . Introduction compensate the effect of the apparent Li enhancement from cool 

spots, when looking at disk integrated spectra. Several studies 
The surface lithium abundance is a crucial diagnostic mark for of the rdation of u abundance with stellar activity were pub- 
studying stellar structure and evolution. The 7 Li isotope, which Hshed xhese studies had different results though9 from a clear 
represents the vast majority of the solar Li, is destroyed through correlation (e . g ., Jeffries et al. 1994, Fernandez & Miranda[[998] 
proton capture at the temperatures of ~3 MK (Lugaro & Chieffi Neuhauser et aL t0 no relation at all, even in the pres- 
mili- This means that the stellar surface layers are progres- ence of very strong starspot activity (e g ? Pa ii avacmi et al. 1993, 
sively depleted of lithium through material transportation to the Mallik 1998 Berdyugina et al 119981 ) 
deeper, hotter layers by convective motions, turbulent diffusion * & 

and meridional flows. Therefore, the surface lithium abundance In a statistical analysis De Medeiros et al. (2000) tested the 
^ , can in principle be used as an age indicator but, in reality, the possible connection between rotation and Li abundance in sin- 
situation is more complex. gle F, G and K-giants. Their result showed a link between the 
Comparison of Li lines observed from the quiet Sun and rotation-rate discontinuity at spectral type G0III and the Li abun- 
sunspot umbra e show s that the L i lines are enhanced in sunspots dance near the same spectral type. Giants cooler than G0III 
(e.g., Grevesse [19681 Giampapa [19841). Thi s is due to the lower are generally very slow rotators (a mean of the vsin/ value of 
temperature of the spots. Giampapa (119841 ) suggested that this 2 kms" 1 ) and generally also show low Li abundances, while 
effect could possibly also be seen in disk-integrated stellar spec- the giants warmer than GO are fast rotators and have higher Li 
tra and that the strength of the lithium line in stars could be abundances. Just recently, Strassmeier et al. (|2012) found a trend 
greatly affected by cool spots. On the other hand, the solar of increased Li abundance with rotation period but with a dis- 
lithium abundance is lower in hot plage regions. This effect could persion of up to 3-4 orders of magnitude. In this study, stars 

in binaries appear to exhibit 0.25 dex less surface lithium than 
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* Based on observations obtained at the Kitt Peak National tation dependent. On the other hand, there is a small group of 

Observatory, U.S.A., and at the Canada-France-Hawaii Telescope, single K-giants which show rapid rotation and also larger than 

U.S.A. expected Li abundances (e.g., Reddy & Lambert 2005 , Drake et 
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al. 120021 Fekel & Balachandran 1993, and references therein). 
Charbonnel & Balachandran (120001 ) concluded that relatively 
high Li abundance is not abnormal for a low-mass star on the 
red-giant branch (RGB). They identified two distinct evolution- 
ary episodes on the RGB in which extra mixing of Li may cause 
an unusual high surface abundance. Both episodes are related to 
the close vicinity of the bottom of the convective envelope to 
the hydrogen-burning shell. Rotational mixing will likely play 
an important role in bringing Li-rich material into the convec- 
tive mantle (cf. Lugaro & Chieffi 2011). Per definition, we call 
a red-giant star "Li-rich" once it has evolved through the first 
dredge-up and still has a logarithmic surface abundance greater 
than 1.5. 

Attention to DP CVn (HD 109703) and DIPsc (HD 217352) 
was drawn in the Vienna- KPNO Doppler-imaging candidate sur- 
vey (Strassmeier et al. 2000). According to this database, the 
stars are rapidly rotating single K-giants, having projected ro- 
tational velocities of ^35 kms -1 (with 2-4 km s" 1 estimated er- 
rors). Charbonnel & Balachandran (120001 pointed out that these 
two stars are new candidates for the small group of rapidly- 
rotating Li-rich giants that appear at or very near the RGB bump 
and already exhibit freshly synthesized lithium. Preliminary de- 
terminations of the Li abundances yielded log^(Li) of 2.82 and 
2.64 for DP CVn and DIPsc, respectively (Strassmeier et al. 
120001) . 

In this paper, we redetermine the absolute dimensions (see 
Sect. 0) and the lithium abundance of these two giant stars 
(Sect. l3.3b . We show them to have indeed a constant radial veloc- 
ity within our precision limits, thus are single stars (Sect. 12. ll) . In 
Sect. 131 we study the surface spot configuration and evolution by 
means of Doppler imaging and complement this in Sect. [5] with 
a long-term photometric spot analysis. Finally, Sect, [^concludes 
on our results. 

2. Observations 

2.1. Spectroscopy and radial velocities 

Most spectroscopic observations were obtained at Kitt Peak 
National Observatory (KPNO) with the 0.9-m coude feed tele- 
scope between March 29 and May 9, 2000 (summarized in 
Table [B. The F3KB 3k-CCD was used together with grat- 
ing A, the long collimator and a 280 /mi slit. This configura- 
tion allows a resolving power of R ^28,000 at a dispersion of 
4.81 kms -1 /pixel. The full wavelength coverage was 300 A and 
the spectra were centered at 6490 A. The exposure time was typ- 
ically 4500 seconds for DP CVn and 1800 seconds for DIPsc, 
giving average signal-to-noise ratios (S/N) of 130:1 and 160:1, 
respectively. 

Six high resolution spectra were obtained with the 3.6m 
Canada-France-Hawaii telescope (CFHT) during May 18 and 
19, 2000. The Gecko spectrograph was used with the 316 1/mm 
grating in 8th order and provided a resolution of 120,000. 
The 4.4kx2k 13.5/mi-pixel EEV1E CCD gives a ^10nm wide 
wavelength range centered at the Li-6708 line. 

Data reductions were done using the NOAO/IRAF soft- 
ware package and followed our standard procedure for coude 
spectra (see, e.g., Weber & Strassmeier 2001 >. Nightly obser- 
vations of various radial-velocity standards were used to ob- 
tain the radial velocities for DP CVn and DIPsc (Table [T]). For 
DP CVn we used mostly 16 Vir (K0.5III, v r =+36.48kms _1 ) as 
the standard except for the spectra at HJDs 2,451,632.65 and 
2,451,633.77 when £Gem (K0III, v r =+3.23kms" 1 ) was used 
and HJDs 2,451,663.75 and 2,451,670.84 when HD 145001 



(G8III, v r =-10.33kms l ) was used. For DIPsc we used mostly 
/?Oph (K2III, v^-mSkms" 1 ) except for HJD 2,451,672.96 
where we used HD 145001, and for the three CFHT spectra 
where we used y Ser (F6IV, v r =+6.58kms _1 ) and (5 CVn (GOV, 
v r =+6.33kms _1 ) (velocity standards are taken from Scarfe et 
al. 119901 ). There appears to be a zero-point shift of 2 kms" 1 
of the DI Psc velocities between KPNO and CFHT data, which 
we interpret to be mostly due to the use of J3 Oph as the KPNO 
template star. Fekel (2004, private communication) obtained five 
more velocities of DI Psc with the coude feed telescope at KPNO 
and got an average velocity of -18.1 kms" 1 , in agreement with 
our CFHT velocities. We therefore assumed the CFHT veloc- 
ities to be correct and accordingly shifted the KPNO data by 
2.0 kms" 1 . The average radial velocity values of DP CVn from 
a total of 23 measurements is -29.6+0.5 (rms)kms" 1 , while 
from 18 measurements the average value for DI Psc is -18.3+0.7 
(rms)kms" 1 . The small rms values suggest that both stars are 
single. 

A summary of the spectroscopic observations is given in 
Tabled] 

2.2. Photometry 

Photometric observations of DP CVn and DI Psc were obtained 
with Wolfgang and Amadeus, the two 0.75-m twin automatic 
photoelectric telescopes (APTs) of the Leibniz-Institute for 
Astrophysics Potsdam (formerly owned by the University of 
Vienna) at Fairborn observatory in Arizona (Strassmeier et al. 
119971) . Amadeus is optimized for red wavelengths and used 
Johnson-Cousins Vic filters while Wolfgang is optimized for 
blue wavelengths and used Stromgren by filters. All the mea- 
surements were done differentially in respect to HD 109529 
for DP CVn and to HD 217019 for DIPsc. Check stars were 
HD 108400 and HD 217590, respectively. For further details 
on the observing procedure and data reduction we refer to 
Strassmeier et al. (119971 ) and Granzer et al. (2001 ). The observa- 
tions of DP CVn consist of a total of 1169 data points, with the 
largest number of data points (381) in the V band. For DI Psc we 
have in total 346 measurements with the largest number (128) in 
the y band. Mean photometric errors by bandpasses are 0™006 
(V), (F007 (7 C ), (T004 (b), and (T003 (y) for both stars in this 
paper. A summary of the photometric observations is given in 
TableEl 

3. Stellar parameters 

3.1. Absolute dimensions 

Hipparcos measured parallaxes for DP CVn and DIPsc 
of 3.35+0.99 mas and 5.27+0.73 mas, respectively (van 
Leeuwen |2007b . This translates into distances of 299 ^j 5 pc for 
DP CVn and 190^3 pc for DIPsc. Table [3] gives a summary of 
the astrophysical data of both systems. 

The brightest observed V and Ic magnitudes for DP CVn are 
8™50 and 7M2 respectively, and for DIPsc 7^10 and 5^95. Due 
to the relatively large distances, interstellar reddening has to be 
taken into account before using absolute magnitudes. The galac- 
tic coordinates of DPCVn and DIPsc are £=128.61 £=+65.45 
and £=78.90 £=-47.96, respectively. We infer the reddening to 
these directions from the dust maps by Schlegel et al. (119981 ) 
and find E(B - V) ^(T014 for DPCVn and E(B - V) ^(F073 
for DIPsc. Using the transformations given by Schlegel et al. 
this translates to A v « 0^05 and A k « 0TO3 for DPCVn, 
and A v « (T24 and A Ic « (F14 for DIPsc. Together with 
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Table 1. Summary of the spectroscopic observations for DP CVn and DIPsc. The table gives the HJDs of the observations, the 
corresponding phases calculated using Eq.[T]and Eq.|2]from Sect. 13.21 respectively, and the measured signal-to-noise (S/N) ratio of 
the spectra, the wavelength region (A/1), the heliocentric radial velocity values with their errors and the observatory. 



DP CVn 


DIPsc 


HJD 


phase 


S/N 


A/i 


Vhelio 


observing 


HJD 


phase 


S/N 


AA 


Vhelio 


observing 


2,451,000+ 






[A] 


[km/s] 


site 


2,451,000+ 






[A] 


[km/s] 


site 


238.92 


0.90 




3895-4110 




KPNO 


068.82 


0.44 




3895^109 




KPNO 


240.85 


0.04 




6470-6787 




KPNO 


070.78 


0.55 




6469-6787 




KPNO 


632.63 


0.00 


145 


6331-6650 


-29.3+0.4 


KPNO 


639.00 


0.00 


133 


6331-6650 




KPNO 


633.75 


0.08 


100 


6331-6650 


-29.6+0.5 


KPNO 


642.01 


0.17 


89 


6331-6650 


-18.6+0.5 


KPNO 


636.66 


0.29 


154 


6331-6650 


-30.4+0.4 


KPNO 


661.00 


0.22 


160 


6331-6650 


-16.5+0.4 


KPNO 


637.75 


0.37 


148 


6331-6650 


-29.3+0.6 


KPNO 


661.99 


0.27 


166 


6331-6650 


-18.5+0.6 


KPNO 


638.70 


0.43 


150 


6331-6650 


-29.7+0.5 


KPNO 


662.99 


0.33 


164 


6331-6650 


-18.0+0.6 


KPNO 


639.79 


0.51 


147 


6331-6650 




KPNO 


663.99 


0.38 


152 


6331-6650 


-18.1+0.5 


KPNO 


640.72 


0.58 


130 


6331-6650 


-29.3+0.4 


KPNO 


665.99 


0.49 


167 


6331-6650 


-17.2+0.4 


KPNO 


641.70 


0.65 


138 


6331-6650 


-29.3+0.4 


KPNO 


666.99 


0.55 


161 


6331-6650 


-18.2+0.4 


KPNO 


642.77 


0.72 


131 


6331-6650 


-29.1+0.5 


KPNO 


667.99 


0.60 


156 


6331-6650 


-18.4+0.6 


KPNO 


644.72 


0.86 


122 


6331-6650 


-28.4+0.7 


KPNO 


668.99 


0.66 


190 


6331-6650 


-18.4+0.7 


KPNO 


660.79 


0.01 


134 


6331-6650 


-29.9+0.5 


KPNO 


669.98 


0.71 


177 


6331-6650 


-18.7+0.6 


KPNO 


661.77 


0.08 


138 


6331-6650 


-29.7+0.5 


KPNO 


670.99 


0.77 


149 


6331-6650 


-17.5+0.7 


KPNO 


662.77 


0.15 


129 


6331-6650 


-29.8+0.4 


KPNO 


671.96 


0.82 


168 


6331-6650 


-18.1+0.5 


KPNO 


663.76 


0.22 


132 


6331-6650 


-29.6+0.5 


KPNO 


672.97 


0.88 


106 


6331-6650 


-19.1+0.7 


KPNO 


664.80 


0.30 


123 


6331-6650 


-29.7+0.4 


KPNO 


673.94 


0.93 


152 


6460-6779 


-18.5+0.6 


KPNO 


665.78 


0.37 


112 


6331-6650 


-29.3+0.6 


KPNO 


673.99 


0.94 


206 


6331-6650 


-18.5+0.5 


KPNO 


666.78 


0.44 


114 


6331-6650 


-29.0+0.5 


KPNO 


683.12 


0.44 


230 


6385-6460 




CFHT 


667.79 


0.51 


118 


6331-6650 


-30.2+0.5 


KPNO 


683.13 


0.44 


230 


6385-6460 




CFHT 


668.79 


0.58 


139 


6331-6650 


-29.6+0.4 


KPNO 


683.13 


0.44 


230 


6385-6460 




CHFT 


669.78 


0.65 


126 


6331-6650 


-29.7+0.4 


KPNO 


684.12 


0.50 


250 


6654-6760 


-19.5+1.7 


CFHT 


670.85 


0.73 


135 


6331-6650 


-29.4+0.5 


KPNO 


684.13 


0.50 


250 


6654-6760 


-18.8+2.1 


CFHT 


671.75 


0.79 


103 


6331-6650 


-29.8+0.5 


KPNO 


684.14 


0.50 


250 


6654-6760 


-18.5+2.2 


CFHT 


672.71 


0.86 


131 


6331-6650 


-30.7+0.9 


KPNO 














673.75 


0.94 


116 


6331-6650 


-30.6+1.0 


KPNO 















Table 2. Photometric observation of DP CVn and DI Psc obtained with the Wolfgang and Amadeus APTs. The table gives the HJDs 
of the observations with the corresponding UT dates and the total number of observations in each bandpass. 



DP CVn 


DIPsc 


HJD 


UT date 


V 


Ic 


b 


y 


V 


Ic 


b 


y 


2451098.7-2451197.6 


12.10.1998-19.01.1999 























89 


2451250.0-2451360.7 


12.03.1999-01.07.1999 











68 














2451434.8-2451539.6 


13.09.1999-27.12.1999 














74 


73 








2451551.9-2451825.8 


08.01.2000-08.10.2000 


104 


97 


134 


133 


16 


14 


38 


39 


2451879.0-2452079.8 


30.11.2000-14.06.2001 


44 


40 


45 


47 














2452177.8-2452179.8 


25.09.2001-27.09.2001 














2 


1 








2452248.0-2452423.7 


04.12.2001-29.05.2002 


80 


81 




















2452629.0-2452786.7 


20.12.2002-27.05.2003 


72 


70 




















2452990.0-2453356.0 


16.12.2003-16.12.2004 


81 


73 






















Total 


381 


361 


179 


248 


92 


88 


38 


128 



the Hipparcos distance this yields an absolute visual magnitude 
My = + r?07+°;^ and the de-reddened V - I c color of 1^06 for 
DP CVn, and M v = +0M7+°;^ and V - / c =r?05 for DIPsc. It 
has to be noted that Schlegei et al.'s dust maps give the inter- 
stellar reddening throughout the entire Milky Way, which means 
that the adopted reddening could be exaggerated in our case. 
Nevertheless, the target stars are located at such distances and 
positions that they are likely above most of the dust, and the 
expected errors are well within the errors caused by the uncer- 
tainties in other parameters. 

The de-reddened V - Ic color can be used to estimate the ef- 
fective temperature of our stars, e.g., with the empirical formula 
given by Mc William (1990) and the transformation of Bessell 
JT979l) . For both stars, we obtain r eff =4600+100K. This is in 
agreement with an independent estimation from our b - y data. 
The bluest observed b — y colors are 0^70 and 0^7 1 for DP CVn 



and DI Psc, respectively, and, with the calibration of Clem et al. 
(12004b , suggest r eff =4550+ 100 K for logg=2.5 for both stars, in 
agreement with the estimate from V - Ic within its uncertainties. 
This gives a bolometric correction of -0.53 (Flower [1996) and 
bolometric absolute brightness of M hol = +(F54+°-^ (DP CVn) 
and M bo i = -(T06^^ (DIPsc). The logarithmic luminosity 
(log ^) of DPCVn is 1.68+^ and of DIPsc is 1.92+°;}j>, based 
on an absolute bolometric magnitude of the Sun of Mb i© = 
4 m 74 (Cox [20001 ) . 

Morel et al. (2004 ) determined temperatures from the excita- 
tion equilibrium of Fe i lines of 13 active binaries and found only 
fair agreement with the temperature from B - V color, but even 
cooler by 200 K on average from V-R and V-I. Despite that the 
temperatures from these three indices agree within their errors in 
case of our two active stars, we employed the line-ratio calibra- 
tion for giants by Strassmeier & Schordan (2000 ). This is nor- 
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DI Psc 



P mt =9.0358d 



5200 5000 4800 
T [K] 



4600 4400 



Fig. 1. The observed positions of DP CVn and DI Psc in the 
Hertzsprung-Russell diagram. The post-main-sequence tracks 
for 1.5, 1.8, 2.0, 2.3, 2.5 and 3.0 solar masses are shown (c.f. 
Granzer et al. 2000). These tracks suggest masses of 1.8 M© 
and 2.3 M© for DP CVn and DIPsc, respectively. Both stars ap- 
pear on the RGB at the end of the first Li dredge-up phase, 
which is marked by the dashed lines taken from Charbonnel & 
Balachandran(20®. 



DP CVn 



Power spectrum 




1400 1600 1800 
HJD (2450000+) 



0.05 0.075 0.1 0.125 



Fig. 2. AV and Ay magnitudes and periodograms from the APT 
data of DP CVn (top) and DIPsc (bottom). The vertical grey 
stripes indicate the time of our Doppler imaging data. The pho- 
tometric periods adopted for this work are 14.010 d (DP CVn) 
and 1 8.066 d (DIPsc, but see also Sect.[X2lfor details). 



mally not done because our Doppler-imaging stars have rather 
strong rotational broadening that makes it virtually impossible 
to measure accurate line-ratios but may serve as an independent 
cross check for the photometric temperatures. Altogether, eight 
line ratios between 6405 A and 6460A were measured for both 
stars and, following the procedure suggested by Strassmeier & 
Schordan (2000, their Fig. 18) for V - I only, convert to av- 
erage effective temperatures of 4530+75 K and 4580+75 K for 
DP CVn and DI Psc, respectively. Note that the DP CVn line ra- 
tios appear 0-3% larger than for DIPsc, which is reflected in 
the different effective temperatures. However, the expected ex- 
ternal uncertainties due to the large line broadening are at least 
of the same order. We conclude that the line-ratio method at least 




DI Psc 



P rot =18.066d 



0.80 
0.82 



0.84 * W * 



S1 



0.80 
% 0.82 
< 0.84 

0.86 



S2 



• • • • • 



0.80 
0.82 
0.84 fc, 
0.86 



S2 




S3 



0.80 
>> 0.82 
< 0.84 

0.86 



S4 



0.80 
0.82 
0.84 : 
0.86 - 



S4 



0.50 
Phase 



0.50 
Phase 



Fig. 3. Folded light curves of DI Psc for four observing runs be- 
tween HJD 2451098.7 and HJD 2451825.8 (S1-S4 from top to 
bottom panels, respectively), assuming rotation period of either 
9.0358 d (left panels) or 18.066 d (right panels), according to 
the power spectrum in Fig. [2] bottom. Note that using 1 8.066 d 
yields less scattered light curves. For a quantitative comparison 
see Sect. 13.21 



agrees with the previous temperature determinations but fails to 
improve the accuracy. 

The positions of DP CVn and DIPsc in the Hertzsprung- 
Russell (H-R) diagram are shown in Fig. [T] We have com- 
puted new post main sequence evolutionary tracks with the 
Kippenhahn code with up-to-date input physics (as described 
in Granzer et al. 120001) . No overshooting was assumed though. 
The comparison with these tracks for solar metallicity suggests 
masses of approximately 1.8+°;J M Q (DP CVn) and 2.3+°;* M Q 
(DIPsc). We placed our stars aiso on the colour-magnitude dia- 
gram constructed with the Padova stellar evolution code (Bertelli 
et al. 2008) (considering overshooting) which yielded slightly 
lower values but still within the given error boxes, i.e., 1 .7+q 3 M© 
for DP CVn and 2.1^ M for DIPsc. This result places both 
stars quite high up on the red-giant branch (RGB) and very close 
to the end of the first dredge-up phase. DP CVn has basically al- 
ready passed this phase. During the first dredge-up the Li abun- 
dance at the surface gets diluted as Li-free material from the 
stellar interior is brought to the surface, i.e., both stars should 
have low Li abundance. 



3.2. Rotational properties 

The rotation period of DP CVn and DI Psc was determined from 
our photometric observations using the frequency analysis pro- 
gram package MuFrAn by Kollath (1990), see Fig. [2[ For both 
of the stars V and y observations were used together, as this 
gave the largest amount of data points for the analysis. The 
most significant period found in the DP CVn observations was 
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Table 3. Summary of astrophysical data. 



Parameter 


DPCVn 


DIPsc 


Distance (pc, Hipparcos) 




190^0 


Spectral type 


K1III 


K1III 


My 


+ l m 07 +0 - 57 

-0.76 


+(T4T° f 2 


Luminosity (log ^) 


1.68^ 


i Q9+0.13' 


log g (adopted) 


2.5 


2.5 


V LQ \ \Ji\^\ U LI LI U 1 1 U LI ) 


1?06 


1?05 


T* (K) 


4600+100 


4600+100 


vsin/ (kms -1 ) 


39.0+1.0 


42.0+1.0 


Rotation period (d) 


14.010+0.040 


18.066+0.088 


Inclination (°) 


35+10 


50+10 


Radius (R Q ) 
Mass (M Q ) 




19.6^ 

9 a +0.4' 
^• J -0.3 


LTE Li abundance (log) 


2.28+0.18 


2.20+0.15 



14.010 + 0.040 days. In the DIPsc data two similarly signifi- 
cant periods were found: 9.0358 + 0.060 d and 18.066 + 0.088 d. 
Although, the shorter period is the more significant one, after 
a careful investigation the longer period was adopted. The rea- 
sons for this are manifold: earlier period analysis (Strassmeier 
et al. E0OO) indicates a period of 18.4d, two other period analy- 
sis methods, the Lomb-method (Press et al. 1992 ) and the Three 
Stage Period Analysis (TSPA, Jetsu & Pelt [T999b also favour the 
longer period. In addition, the individual light curves (S1-S4) 
in Fig. [3] show less scatter when phasing them with the longer 
period. For a quantitative test we apply spot model fits for the 
seasonal light curves (SpotModeL by Ribarik et al. 120031) to 
compare the sum of squares of residuals, assuming either P rot 
of 9.0358 d or 1 8.066 d. Unambiguously, systematic improve- 
ment of the fits is found when using the longer period, resulting 
in 53% (SI), 30% (S2), 89% (S3) and 12% (S4) less residu- 
als. Finally, we note that in a situation where there are two spot 
groups separated by approximately half a rotation, many period 
search methods would prefer the period which is half of the real 
one. 

Rotational phases were calculated using the following 
ephemerides: 

HJD = 2,451,632.63084+ 14.010(40) x E, (DPCVn), (1) 

HJD = 2,451,639.00330+ 18.066(88) xE, (DIPsc) (2) 

where the epoch is the time of the first spectroscopic observa- 
tion in the data set and the period is the photometric period. The 
errors of the period determinations are estimated by increasing 
the residual scatter of the nonlinear least- squares solutions a cer- 
tain degree, which corresponds to 10% of the accuracy of our 
photometric data (cf. Olah et al. 120031) . 

Our Doppler imaging analysis (Sect. HI) helped us to deter- 
mine the projected rotational velocities (vsin/) and the incli- 
nation angles i from goodness of fit landscapes (see Fig. [6]). 
When we combine the resulting vsin / of 39.0+1.0 kms -1 for 
DPCVn and 42.0+1.0 kms" 1 for DIPsc with above rotational 
periods, the minimum radii of our stars are 10.7+0.3/?© and 
15.0+0.47?©, respectively. Doppler imaging constraints the in- 
clination to 35°+10° for DPCVn and 50°+10° for DIPsc, giv- 
ing absolute radii of 18.77?© and 19.67?©, respectively (see also 
Table 0). We note, that within the given errors, the radii and 
masses from Table [3] are in agreement with the adopted log g 
value. 
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Fig. 4. Top panel: Lii 6708 A spectra of DIPsc (top) and 
DPCVn (shifted by -0.2). Also shown is a spectrum of the in- 
active K0.5III giant 16 Vir (shifted to bottom) for comparison 
purpose. The measured equivalent width of the lithium line of 
DIPsc and DPCVn is 328+5 mA and 316+7 mA, respectively. 
Bottom panels: Lii 6708 A line fitted by synthetic spectra for 
DPCVn (left) and DIPsc (right), resulting in logarithmic abun- 
dances of log?r(Li)=2.28, and of 2.20, respectively. 



3.3. Lithium abundances 

Top panel of Fig. |4] shows a high-resolution spectrum of the ex- 
ceptionally strong Li i 6708A line for both stars in comparison 
with the inactive K0.5III star 16 Vir (otherwise an IAU velocity 
standard; see, e.g., Strassmeier et al. 1998 Weber & Strassmeier 
2001, Strassmeier et al. 120111 etc.). The measured equivalent 
widths are 316+7 mA (DPCVn) and 328+5 mA (DIPsc). The 
corresponding equivalent width (EW) of the sum of the lines 
that would appear as blends within the Li line width, as mea- 
sured in 16 Vir, is 39.8+3 mA. The uncertainty of just 5-7 mA 
is estimated from different types of fits to the line profile, e.g. 
with a Gaussian profile of variable width where only the profile 
points of the red profile wing are used for the least- squares fit, 
from a normal Gaussian fit or from a simple profile-area integra- 
tion. The final Li EWs for the two target stars are obtained by 
subtracting the 16 Vir contribution and are 277mA for DPCVn 
and 288 mA for DIPsc. 

Next we fit synthetic (LTE) spectra to the Li line profiles. 
Synthetic spectra were calculated and matched to the observed 
spectra by adjusting the Li abundance to fit the 6707.8 A feature. 
We used the VALD database (Piskunov et al. |19951 Kupka et al. 
119991) and Kurucz's CD-ROM No. 23 for the numerous atomic 
and molecular lines. Its log gf values were slightly tuned using 
a spectrum of f3 Gem, while we adopted (and fixed) the values 
from Reddy et al. (2002 ) with the component structure of the Li 
resonance doublet from Hobbs et al. ( 1999 ). We assume that Li 
is purely 7 Li. 

Bottom panels of Fig. |4] show comparisons of the observed 
and synthetic spectra for both stars. Li abundances derived from 
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Fig. 5. Ua line profiles of DP CVn (left) and of DIPsc (right). 
In the top panels the observations are plotted together, while in 
the bottom panels the corresponding averages are subtracted to 
emphasize the nightly variations. Note the telluric water lines, 
especially at around 6552.5 A, 6557 A, and 6564 A, as well as 
the sporadic cosmic peaks. 



our LTE analysis are equal to 2.28 for DP CVn and 2.20 for 
DIPsc. Small perturbations of the observed profiles are due to 
spots on the stellar surfaces which may influence the abundance 
value. Note that 100 K uncertainty in T e ff yields ^0.2 dex. 

Tests with changes of the continuum level of 2% resulted 
in abundance changes of 0.10-0. 12 dex. In our case the contin- 
uum level of the observed spectra was obtained iteratively from a 
comparison with calculated spectra including rotationally broad- 
ening. Notice a 4% depression of the continuum in the bottom 
panels of Fig.|4ldue to the rotational convolution with several ab- 
sorption blends. We did not carry out a full abundance analysis 
and thus assumed solar abundances for all other elements except 
Li and C,N,0. As our target stars likely belong to the "bump" 
clan (Charbonnel & Balachandran 2000) corresponding values 
of C,N,0 abundances were used in our calculations. 

Non-LTE corrections interpolated from the results by 
Carlsson et al. (119941 ) were found to be small (<0.03dex) 
and comparable to its uncertainties and methods of inter- 
polation. Therefore, we refrain from giving non-LTE abun- 
dances. However, we note that the non-LTE curves of growth 
of Pavlenko & Magazzu (fT996b for a 4600K/logg=2.5 model 
convert the EWs measured into a Li abundance of 2.37+0.05 
(DP CVn) and 2.40+0.05 (DIPsc), both on the log n(H)= 12.00 
scale, i.e., larger by 0.1-0.2 dex than our LTE values. 

The position in the H-R diagram in Fig. Q] and the large Li 
abundance places both stars at the end of the first dredge up ac- 
cording to the models of Charbonnel & Balachandran (2000). 
Thus, they are truly Li-rich stars. They even clump together 
with stars having freshly synthesized Li on their surface, such 
as the super-meteoritic Li stars HD 9746 and HD 233517 (Fekel 
& Balachandran[l993l). 

3.4. Ha line profiles 

In both stars Ha appears as a relatively stable absorption line 
without periodic variations from one spectrum to the next 
(see Fig. 0, suggesting a rather homogeneous chromosphere 
made up of symmetrically distributed emission regions, com- 
parable to a more active version of the solar chromospheric 
network. The averaged EWs are 1.183 A for DP CVn and 
1.369 A for DIPsc, virtually without any significant phase 



dependency. The only markable phase independent change 
is seen in three spectra of DP CVn taken in three subse- 
quent nights from HJD 2,451,636.657 (0=0.288) through HJD 
2,451,637.745 (0=0.365, with EW of 1.526 at maximum) to 
HJD 2,451,638.698 (0=0.433). The line is broadened asymmet- 
rically towards the blue wing, however, the feature did not re- 
turn during the next rotation cycle. Such a non-recurring phe- 
nomenon can be caused by an erupting coronal prominence. A 
similar change in the Ha line was reported e.g., for the effec- 
tively single G8-giant CM Cam by Strassmeier et al. ( 1998 ). 

4. Doppler imaging 

4.1. The TempMap inversion code 

A full description of the TempMap code can be found in Rice et 
al. ([T9891 Rice & Strassmeier (l2000t and in Rice (12002b . The 
program performs a full LTE spectrum synthesis by solving the 
equation of transfer through a set of ATLAS-9 (Kurucz 119931 ) 
model atmospheres at all aspect angles and for a given set of 
chemical abundances. Atomic line parameters are obtained from 
VALD (Piskunov et al. fT995l Kupka et al. [T999l) . Spectra were 
calculated for temperatures ranging from 3500 K to 5500 K in 
steps of 250 K and with solar abundances. Simultaneous inver- 
sions of the usual mapping lines as well as of the photometric 
light variations in two photometric bandpasses are then carried 
out using maximum-entropy regularization. 

A large number of trial inversions helped us to fine tune 
some input parameters, and led us to the final values for the 
projected rotational velocity v eq sin i and the inclination i, which 
were determined with test runs by minimizing the x 2 > i.e., the 
goodness-of-fit of the Doppler reconstructions (for the x 2 defi- 
nition see Rice & Strassmeier 2000 ). The goodness-of-fit land- 
scape in Fig. [6] demonstrates the change of x 2 over a meaning- 
ful domain of the v eq sin i - i parameter plane for both stars. 
Best values for v eq sin/ were 39. 0+ 1.0 km s" 1 (DP CVn) and 
42.0+l.Okms -1 (DIPsc), and for the inclination were 35°+10° 
(DP CVn) and 40°+ 10° (DIPsc). The best combinations are 
listed also in Tabled 

4.2. Doppler reconstructions for DP CVn 

Data for this star were taken in the modus of 12 nights on, 16 
nights off, 13 nights on, i.e., approximately two rotation cycles 
were covered by the observations, with one rotation between 
them. This enabled two independent data sets with one spectrum 
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Fig. 6. The dependency of the normalized goodness of fit (x 2 ) of 
our Doppler maps over the v sin i-i parameter plane. Best values 
forDPCVn (left) are v sin /of 39.0+l.Okms" 1 with/ * 35°+10° 
and for DIPsc (right) are 42.0+ 1.0 km s" 1 with 50°+ 10°, respec- 
tively. 
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Fig. 7. Surface temperature distribution maps of DP CVn for the two available data subsets taken in April 2000 (serl, top row) 
and in May 2000 (ser2, bottom row), using Cai 6439 A, Fei 6430 A, Fei 6421 A, and Fei 6400 A mapping lines. The temperature 
maps are presented in a pseudo-Mercator projection from latitude -35° to +90°. The phases of the observations are marked by 
arrows underneath. Below surface temperature reconstructions, in the middle panels the fitted line profiles are shown. Small dashes 
represent the data points and measure the +l-cr errors. Simultaneous light curves in Stromgren b and y are plotted in the bottom 
panels with their respective fits. 



per night and a total of 10 spectra for the first rotation (serl, 
note that one Cai 6439 A spectrum was omitted due to probable 
cosmic noise) and 14 for the second rotation (ser2). Four nights 
were lost during the first rotation due to bad weather, the second 
rotation was covered completely. 

The sampling of our spectra allowed ^23 pixels across the 
full width of a given spectral line and the S/N per pixel is be- 
tween 100:1 to 150:1 as measured with IRAK Additionally, 
there are 17 simultaneous two-colour photometric data points 
in Stromgren b and y. 

Maps for the two data subsets taken in Apr 2000 (serl) and 
in May 2000 (ser2) are shown in Fig. [7] for Cai 6439 A, Fei 
6430 A, Fei 6421 A, and Fei 6400 A, our most commonly used 
absorption lines with well-known formation physics and mostly 
free of telluric blends. (Surface reconstruction from Li can not 
be done because we have only one high-resolution Li spectrum.) 
Also shown are the simultaneously fitted light curves in b and y. 
The individual reconstructions agree quite well, revealing dom- 
inant spots at lower latitudes with ^700 K below T e ff, as well 
as some high-latitude features but mostly with lesser contrast. 
Note also, that the minimum temperature value of the Ca-maps 
in Fig. [7] is set to 4200 K, i.e., 300 K hotter than that of the Fe- 



maps. The less contrast of the Ca-line reconstruction is proba- 
bly due to the fact that the Ca line has broader equivalent width 
and also is more sensitive to temperature changes than the iron 
lines. On the other hand, smaller differences could occur through 
the relatively low S/N of the spectra. The agreement of the indi- 
vidual reconstructions can be demonstrated by comparing them 
with their averages (see Fig. [8}. Markedly, most of the notice- 
able features in the contemporaneous individual maps in Fig. [7] 
can be identified in the corresponding combined map. 

4.2.1 . Searching for surface differential rotation 

Latitude dependency of the spot displacements from the first 
subset (serl) in Apr 2000 to the second one (ser2) in May 2000 
can be assigned as a trace of surface differential rotation. Such a 
dependency is suspected when looking at Fig. [TJ e.g., where the 
spots near the equator seem to be shifted mostly backwards in 
longitude, unlike spots at higher (>45°) latitudes. Differential ro- 
tation can be measured by cross-correlating subsequent Doppler 
images (Donati & Collier Cameron 1 1997I ). We apply the cross- 
correlation method called ACCORD, as described recently in 
Kovari et al. (12012b . We cross-correlate the corresponding lon- 
gitude strips of the subsequent Doppler maps for each latitude 
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value (in practice in bins of 5° width) for the four available indi- 
vidual reconstructions. This way, we get four cross-correlation 
function (CCF) maps. Then, the CCF maps are combined to 
achieve the average correlation pattern, wherein the sign of the 
differential rotation is expected to be intensified as a common 
feature of the individual CCF maps. Finally, the best correlated 
pattern is fitted by a solar-like differential rotation law, assuming 
the usual quadratic form of Q(J3) = Q e q _ sin 2 /?, where Q e q 
is the angular velocity at the equator and AQ = Q eq - Q po i e is 
the difference between the equatorial and the polar angular ve- 
locity. The dimensionless surface shear parameter is expressed 
as a = AD/Deq- 
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Fig. 8. Combined images of DP CVn for the two data sub- 
sets serl (top) and ser2 (bottom). The maps are presented in 
a pseudo-Mercator projection from latitude -35° to +90°. The 
phases of the observations are marked by arrows below the maps. 
Note, that even if the individual reconstructions show differences 
(cf. Fig. [7]>, the average maps still retain most of the traits. 

In Fig.[9]the best-fit rotation law to the correlation peaks per 
latitude bin (represented by fitted Gaussian maxima and their er- 
rorbars) is plotted. The relevant part of the correlation pattern 
extends from -25° to +40°, and is strongest at around 30°, i.e., 
the latitude belt that rotates with the disk-integrated photomet- 
ric period. This is not surprising because below the equator the 
inclination angle limits the visibility while at higher latitudes 
the lack of strong surface features that are jointly present on the 
corresponding maps yields weak and slurred correlation. Note 
that at even higher latitudes, close to the visible pole, the ge- 
ometric singularity makes "correlation" comparably meaning- 
less. The resulting fit indicates antisolar-type differential rota- 
tion, i.e. equatorial deceleration, with £2 eq = 25.336 ±0.044°/day 
and AQ = 0.875 ± 0.411°/day or a = -0.035 ± 0.016 (errors 
are estimated from the FWHMs and amplitudes of the Gaussian 
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Fig. 9. Cross-correlation function (CCF) map for DP CVn. The 
image is the average of four CCF maps obtained from cross- 
correlating the individual Doppler reconstructions (i.e., Cai 
6439 A, Fe 1 6430 A, Fe 1 642 1 A, and Fe 1 6400 A maps) taken in 
Apr (serl) and May (ser2). Dots with error bars are the Gaussian- 
fitted correlation peaks per 5° latitude bin. The best fit suggests 
an anti-solar differential rotation function (solid line) and corre- 
sponds to a surface shear a of -0.035. 



profiles). We note, that considering the moderate significance, 
the surface shear parameter needs further validation from better 
quality data. 



4.3. Doppler reconstructions for D I Psc 

Our observations allow for one Doppler reconstruction from data 
taken between HJD 2,451,661.00 and HJD 2,451,673.99, i.e. 
covering 13 nights, or 72% of one rotation cycle. Unfortunately, 
no simultaneous photometric data were available (but see 
Strassmeier et al. |1993l on the role of photometry in TempMap). 
Fig. [TO] shows the Cai 6439 A, Fei 6430 A, and Fei 6400 A im- 
ages along with the corresponding spectroscopic data and their 
respective fits. 

Attempts to use other mapping lines from the observed spec- 
tral region, i.e. Fe i 6421 A, Fe i 641 1 A or Fe 1 6393 A led to un- 
acceptable results probably due to unknown blends (atomic and 
molecular), to uncertain turbulence parameters, etc. Note that the 
available Li spectra of DI Psc are insufficient for Doppler recon- 
struction. 

The Cai 6439 A map reveals cool spots at high latitudes 
as well as spots at lower latitudes, but with lesser contrast. 
Similarly, the Fe 1 6430 A, and Fe 1 6400 A maps show cool polar 
features, but dominant spots appear also at latitudes below ^40°, 
forming a belt-like structure. The temperature scale is different, 
typically ranging from AT" « 800 K in the Ca map up to 1300 K 
in the iron maps. Aside from that and from minor dissimilarities 
(most likely due to the moderate quality of the data and the un- 
covered phase gap from ^340° to ^80°), the agreement is still 
acceptable, as shown in Fig. [TT] where an average map is com- 
bined from the individual reconstructions. It is clearly worth to 
revisit this star with instruments with larger wavelength cover- 
age and further explore its surface inhomogeneity. 
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Fig. 10. Reconstructed Doppler maps of DI Psc for Ca 1 6439 A, Fe 1 6430 A, and Fe 1 6400 A mapping lines. The maps are presented 
in a pseudo-Mercator projection from latitude -50° to +90°. The phases of the observations are marked by arrows underneath. 
Below surface temperature reconstructions the fitted line profiles are shown, where small dashes represent the data points and 
measure the + l-cr errors. 



5. Long-term photometry of DP CVn 

5.1. Spot evolution 

The spot evolution was investigated first by plotting the maxi- 
mum, minimum and mean magnitudes in V and y at different 
observing periods against the amplitude. No systematic trend in 
the mean magnitude with amplitude is seen, while the maximum 
and minimum magnitudes are systematically higher and lower, 
respectively, when the amplitude increases. We conclude that on 
DP CVn the light curve variations are caused by a rather sym- 
metrical rearrangement of a more or less constant amount of dark 
spots. This is consonant with the results from Doppler imaging 
(Sect. I4.2K where the most dominant spots are located within a 
belt below ^40° and distributed more or less evenly (see Fig. [8). 
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Fig. 11. Combined image of DIPsc from the available three in- 
dividual Doppler reconstructions. The map is presented in a 
pseudo-Mercator projection from latitude -50° to +90°. The 
phases of the observations are marked by arrows below the map. 



5.2. Active longitudes 

The photometric dataset of DP CVn is suitable for carrying out a 
more detailed analysis of the spot evolution through the 6 year- 
long observing range. For this purpose the V and y observations 
were divided into 16 subsets. These individual light-curves were 
then used for obtaining spot filling factor maps of the stellar sur- 
face. For the inversions we use the Occamian approach method 
(B erdvu gina 1 1 99 8 1 but see also Lanza et al. 119981 ). The stellar 
surface was divided into areas of 10° x 10°. The temperatures 
for the unspotted surface and the spots were fixed to 4600 K and 
4000 K, respectively, based on the Doppler maps (see Sect. 14. 21) . 
The limb-darkening coefficient for the V band for the unspotted 
surface was estimated to be 0.802 from the results by Al-Naimiy 
(1978 ) and the same value was also adopted for the spotted areas. 
As unspotted magnitude we used 8™45 as an estimation from the 
brightest observed magnitude of 8™50. 

The inversions of the individual light curves covering the in- 
terval 1999-2005 are presented in Fig. [l2]together with the light 
curves themselves. As can be seen, the maps usually show two 
active regions which often are extended or even merged. Also, 
these active regions seem to be quite stably separated by half 
of the phase. This kind of activity pattern has been found on 
many active stars (see, e.g., Olah et al. 1991, Henry et al. 1995, 
Jetsu |1996l Berdyugina & Tuominen |1998[ etc.), and may bear 
also some evidence for the existence of a flip-flop phenomenon, 
where the dominant part of the spot activity shifts 180° in longi- 
tude over a short period of time and remains at this new active 
longitude for some time, usually a few years (see e.g., Jetsu et al 
[T9951 Berdyugina & Tuominen [T998l Korhonen et al. l200Tl) . 



5.3. Differential rotation 

Differential rotation can be assessed from seasonal photomet- 
ric period changes, since the current period depends on the 
stellar latitude of the spots that produce rotational modula- 
tion. Therefore, we derived rotation periods for the ten avail- 
able observing seasons (cf. Fig. |2]). The values fall between 
P™ n =13.61d (2nd season) and />™ x = 14.39 d (7th season). A 
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Fig. 12. Light curve inversions for DP CVn. Darker regions indicate larger spot filling factors. The grid in the maps indicates the 
equator and four longitudes separated by 90°. The observations (crosses) and calculated V and y light curves (solid lines) are 
also presented. Note that for convenience reasons we use 2,451,250.99805 as the zero phase in the plot, i. e., th e time of the first 
photometric observation we have, yielding a phase shift of +0.240 as compared to the Doppler maps in Sect. 14.21 ). 



rough guess of the surface shear parameter a can be estimated 
from AP/P, getting ^0,06. We underline nonetheless, that the 
sign of the surface shear (i.e., solar- type equatorial acceleration, 
or opposite) cannot be determined without additional informa- 
tion. This result is quite close to the value derived from spec- 
troscopy. On the other hand, the typical error of the period deter- 
mination of 0.36 d for such short intervals makes this determina- 
tion ambiguous. 

To perform a more reliable analysis we divided the whole 
photometric dataset into two parts. During the first half of the 
observing range (lst-5th seasons) the seasonal amplitudes were 
smaller, while the mean brightness level was higher, as com- 
pared to the second half (6th- 10th seasons) where the ampli- 
tudes were increased and the mean brightness was decreased. 
This change can be interpreted by shifting dominant spotted re- 
gion(s) towards lower latitudes. If so, we would expect longer 
rotation period for the second half of the dataset when assum- 
ing antisolar-type differential rotation. Fourier analysis resulted 
inP l m = 13.862±0.035dforthe first term (lst-5th seasons), and 
Pf ot = 14.012 ± 0.028 d for the second term (6th- 10th seasons), 
implying a minimum shear of a = -0.01. We conclude that 
such an overall period change is in agreement with the proposed 
antisolar-type differential rotation law for DP CVn in Sect. 14.2.11 

6. Conclusions 

The revised Hipparcos parallax places both DP CVn and DI Psc 
at a position in the H-R diagram where a deep convective enve- 
lope coexists with a hydrogen-burning shell. In that evolutionary 
phase, rapid redistribution of the moment of inertia must occur 
or had occurred (e.g. do Nascimento et al. 2000 ) and is likely to 
be a time of high dynamo activity due to the extra shear forces. 
Indeed, a number of overactive single giants were identified from 
their unusual high Li abundance (e.g. Brown et al |1989l Fekel & 
Balachandran l 19931) . Therefore, these stars seem to have already 
experienced their first Li dredge-up according to the models of 
Charbonnel & Balachandran (2000). Both stars of the present 
study likely have developed a degenerate helium core, and will 
later undergo a He-flash, and will evolve on to the early AGB 
and experience a second Li dredge-up. 

The surface spot distribution in both stars is rather unspec- 
tacular in that there is no large polar cap, or a single high-latitude 
feature or any other dominating feature. This is qualitatively in 
agreement with the absence of rota tiona l modulation in Ha at 
that time (for both stars, see Sect. I3.4K and suggestive of an 
atmosphere probably entirely penetrated and saturated by mag- 



netic fields. However, this may not be the case at all times, as 
shown by the light curve inversion for DP CVn in 2002.26, when 
spots were concentrated just in one hemisphere. 

Both stars are very similar to the single high-Li K2-giant 
HD 31993 (Strassmeier et al. 2003 ) in that they are very active 
and in an evolutionary episode where we might expect a rela- 
tion between its magnetic surface activity and the Li abundance. 
Knowing the surface Li distribution and its relation to chro- 
mospheric activity indicators would be most useful for abun- 
dance determinations. Our Doppler maps of HD 31993 and the 
two stars in this paper may support the analogy with the Li en- 
hancement in sunspots, i.e., that the high Li abundance of these 
stars is related to the increased and nearly isotropic spottedness. 
Unfortunately, the quality of our spectra is only moderate and 
higher S/N ratio and higher spectral resolution are warranted. 
Also, further Doppler images of Li-rich stars are needed to quan- 
titatively elaborate on such a speculation. 
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